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We argue that in at least a portion of the history of the universe the relic background neutrinos 
are spatially-extended, coherent superpositions of mass states. We show that an appropriate quan- 
tum mechanical treatment affects the neutrino mass values derived from cosmological data. The 
coherence scale of these neutrino flavor wavepackets can be an appreciable fraction of the causal 
horizon size, raising the possibility of spacetime curvature-induced decoherence. 
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In this letter we point out a curious feature of the neu- 
trino background and the potential implications of this 
both for cosmology and for neutrino physics. Neutri- 
nos and antineutrinos in the early universe should be 
in thermal and chemical equilibrium with the photon- 
and e^-plasma for temperatures T > T^oak ^ 1 MeV. 
For T <ti Tweak, the neutrinos and antineutrinos will be 
completely decoupled, comprising seas of "relic" particles 
freely falling through spacetime with energy-momentum 
and flavor distributions reflecting pre-decoupling equilib- 
rium plus the expansion of the universe. This is analo- 
gous to the presently-decoupled cosmic microwave back- 
ground (CMB) photons. These photons had been cou- 
pled and in thermal equilibrium in the early universe 
when T>Tem~ 0.26 eV. 

Experiments have demonstrated that the neutrino en- 
ergy (mass) eigenstates {vi) are not coincident with the 
weak interaction (flavor) eigenstates {va)- These bases 
are related by the Maki-Nakagawa-Sakata (MNS) matrix, 
|i^a) = J2t^aiWi)^ where a = e,fi,T, and i = 1,2,3 de- 
notes the mass eigenstates, with corresponding vacuum 
mass eigenvalues m^, and Uai are the unitary transforma- 
tion matrix elements (parameterized by 3 mixing angles 
(^12, (^23, (^13, and a CP-violating phase 6). 

Absolute neutrino masses remain unknown, but the 
neutrino mass-squared differences are measured. Stud- 
ies of atmospheric neutrinos reveal a characteristic mass- 
squared splitting, |(5m^tj„| « 2.4 x 10"'^ eV^, associated 
with ^ Vr mixing with sin^ 623 ~ 0.50. Likewise, so- 
lar neutrino observations and reactor-based experimen- 
tal data suggest that the solar neutrino deficit prob- 
lem is solved by flavor transformation in the i>e ^ ^fj-fr 
channel and that the characteristic mass-squared differ- 
ence for this solution is SitLq w 7.6 x 10~^ eV^, and that 
sin^ 012 ~ 0.31. Current experimental limits on ^13 are 
sin^6'i3 < 0.040 (2cr limit) []]. The CP-violating phase 
and ordering of the neutrino mass splittings remain un- 
constrained by experiment. In the "normal" neutrino 
mass hierarchy the solar neutrino mass-squared doublet 
lies below the atmospheric doublet; in the "inverted" 
mass hierarchy it is the other way around. This is il- 
lustrated in Fig. [1] 

For T > Tweak, prior to decoupling, solutions of the 
quantum kinetic equations for neutrino flavor evolution 
show that neutrinos will be forced by weak interaction- 
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FIG. 1: Neutrino mass eigenvalue ordering in the "normal 
mass hierarchy" at left; "inverted mass hierarchy" at right. 



mediated scattering into flavor eigenstates 0, H, 01 ■ At 
epoch Tweak, Corresponding to scale factor Cweak, the neu- 
trino distribution functions for each flavor will be two- 
parameter, Fermi-Dirac black bodies, with the number 
density of Va 's in energy interval dE^ given by 
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where we employ natural units with Ti = c = ks — 1, 
the degeneracy parameter (ratio of chemical potential to 
temperature) for neutrino species Va is rjij^ , and where 
the neutrino energy-momentum dispersion relation is 



,1/2 



(2) 



with p = IpI the magnitude of the spacelike momentum. 
Here we will ignore small spectral distortions between Vf, 
and Ufj^ r from e^-annihilation. 

Subsequent to decoupling, these pure flavor states can 
be regarded as coUisionless, simply free falling through 
spacetime with their momentum components redshifting 
with the scale factor, p oc . As a result, the number 
density of Uq, at an epoch with scale factor a is 
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where T^(a) — TwoakOweak/a is an effective "tempera- 
ture" of the neutrinos, e — p/Ty{a) is a co-moving invari- 
ant, and the energy-momentum dispersion relation is 
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In both of these expressions, we have neglected correc- 
tions of order (mi/Tweake)^ which are small since ex- 
perimental measurements of neutrino masses have ^ 
1 eV, Tweak ~ 1 MeV, and small values of e are sup- 
pressed in the neutrino distribution function, Eq. Q (i.e., 
(m,;/Twoake) ^5 10^^''). The energy-momentum disper- 
sion relation, Eq. ([4]), leads to the standard choice for the 
effective mass of a neutrino in flavor state z/q, 

i 

This is the dynamical mass for an ultrarelativistic neu- 
trino of flavor i>a- However, we will show that when 
the neutrino momentum redshifts to a point where the 
neutrino's kinematics become less relativistic 5], this ef- 
fective mass is no longer relevant in characterizing the 
energy-momentum dispersion relation. 

Examining Eqs. ([3|) and ([4]), we notice that the dis- 
tribution function for Va maintains a self-similar form 
with e and rj^^ being co-moving invariants. This form is 
suggestive of, though subtly different from a Fermi-Dirac 
blackbody with temperature T^, and degeneracy param- 
eter r]^^ . The energy distribution functions for neutri- 
nos in mass eigenstates are a weighted sum of the flavor 
eigenstates, 

dn^^ ='^\Uai\'^dTi^^. (6) 

a 

If the degeneracy parameters for the three flavors were all 
equal, then the distribution function for the mass states 
would have the same Fermi-Dirac form. This follows from 
unitarity, J2a Wai\'^ = 1- However, in general the distri- 
bution functions of neutrinos in mass eigenstates would 
not have a Fermi-Dirac form where the degeneracy pa- 
rameters were not identical for all three active flavors. 

CMB- and large scale structure-derived neutrino mass 
limits typically are predicated on an assumption that 
these distribution functions have a Fermi-Dirac black 
body form. This assumption is invalid when the energy 
spectra of the various neutrino flavors are not identical, 
i.e., when the neutrino degeneracy parameters (or, equiv- 
alently, the corresponding lepton numbers) are unequal. 

The measured solar neutrino mixing parameters cou- 
pled with Big Bang Nucleosynthesis (BBN) considera- 
tions limit r]^^ < 0.15 [1, [7| for all flavors and dictate 
that the degeneracy parameters be within an order of 
magnitude of each other d, 0, [1]. Consequently, the 
spectral and number density differences for the neutrino 
flavors stemming from disparate lepton numbers cannot 
be large, but conceivably could be < 15%. Nevertheless, 
as neutrino mass limits improve with future higher pre- 
cision CMB experiments and better large scale structure 
observations, it may be necessary to include the possibil- 
ity of non-Fermi-Dirac energy spectra in these neutrino 
mass analyses. For the remainder of this letter, we will 
assume that all the neutrino degeneracy parameters are 
zero to illustrate another important result. 
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FIG. 2: [pi Pv)/ pv as a function of redshift. Solid curve is 
for mi — 1 meV in the normal neutrino mass hierarchy and 
the dashed curve is for ms = 1 meV in the inverted hierarchy. 

There is a more general issue associated with deter- 
mining neutrino masses from cosmological data. A naive 
method of calculating the energy density of neutri- 
nos in the universe, is to use the effective mass in Eq. 
([5]) and the number density distributions of neutrinos in 
flavor eigenstates, 

pi"^=E /(^''+-cff,.J^/'rf-..- (7) 

However, if we wish to calculate the energy density of 
these quantum mechanical particles, we ought to consider 
the energy (mass) eigenstates in the calculation, 

i 

= f |C^«P(p' + m2)i/2rf„^^. (8) 

At late times (large scale factors or equivalently small 
redshifts), the momenta of the neutrinos have redshifted 
to a point where they are comparable to the neutrino 
masses. In this case the quantum mechanically consis- 
tent calculation diverges from the naive effective mass 
approach. Fig. [2] shows the fractional difference between 
these two calculation techniques. At large redshifts, neu- 
trino momenta are large enough that all three mass eigen- 
states are ultrarelativstic and masses have little effect on 
the overall energy density. At low redshifts, the kine- 
matics of the neutrinos is no longer ultrarelativistic and 
rest mass becomes significant in the overall energy den- 
sity. At this point the difference in the two methods 
manifests itself in the difference between the quantum 
mechanical expectation value for the mass, (m), and the 
root-mean-squared mass expectation value, {rr?)^^'^ . 

As can be seen in Fig. [21 the disparity in calculating 
the value of becomes significant, of order 10%, for 
redshifts of order 1 — 10, at least for nii ^ 1 meV. This 
could be important because this is the epoch of struc- 
ture formation in the universe where the energy density 
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FIG. 3: The lightest mass eigenvalue as a function of . 
Solid curves use the quantum mechanical calculation method 
described here, while dashed curves use the naive method 
illustrated in Eq. ([7]). Curves on the left correspond to the 
normal neutrino mass hierarchy; curves on the right are for 
the inverted hierarchy. Also plotted is the limit from the 
WMAP 5-year data, Q,^h^ < 0.0071 (95% confidence limit). 



in neutrinos and the character of their kinematics are 
relevant to the formation of large scale structures. 

One method of experimentally determining the mass of 
the neutrino using the CMB is by inferring the transfer 
function in the matter power spectrum at large wavenum- 
bers (or small scales). Massive neutrinos contribute to 
the closure fraction in cold dark matter at the current 
epoch (as inferred by measurements of the CMB power 
spectrum), but at higher redshifts these neutrinos may 
act as hot dark matter. As a result, the amount of cold 
dark matter at these earlier epochs is reduced, suppress- 
ing the formation of large scale structure compared to a 
situation without massive neutrinos. A comparison of the 
observationally-inferred matter power spectrum with the 
power spectrum expected without the effects of massive 
neutrinos implies fi^, the closure fraction contributed by 
neutrinos. Determining the energy density in neutrinos 
at late times is an important aspect of this procedure and 
is necessary in deriving neutrino masses. 

This is a concern since the naive computation of neu- 
trino energy densities differs from the quantum mechan- 
ical method. Fig. [3] shows the relationship between the 
lightest mass eigenvalue (mi in the normal hierarchy or 
ma in the inverted hierarchy) and the measured value 
of We see that the limit on ^^h? for excluding the 
inverted hierarchy is reduced by about 15%. (Here h is 
the Hubble parameter in units of 100 km s"'^ Mpc~^.) 
In addition, at low mass values, we see that the mass 
is a very sensitive function of Q,^h? . As a result, mass 
measurements using this method are sensitive to the pre- 
cision in the inferred measurement of Q,^h? . A 10% error 
in the measurement of Q,^h? would prevent us from in- 
ferring a lightest neutrino mass ^ 5 meV. For neutrino 
masses ^ 10 — 20 meV, there can be a significant error in 
determining the lightest mass eigenvalue if the quantum 
mechanical treatment we present here was not used. 



The current observational limits (shown in Fig.[3]as the 
five-year WMAP data |9|]), are nowhere near the regime 
where a consistent quantum mechanical treatment is nec- 
essary. However, the quantum mechanical analysis may 
be relevant for higher precision future neutrino mass in- 
ferences based on transfer function arguments. 

One possible complication has not been discussed. It 
has been claimed that the phase space density of the relic 
neutrinos is high so that quantum coherence between 
momentum states may develop as a result of fermion 
exchange symmetry p^]. We note that the initial dis- 
tribution of neutrinos, Eq. ([1]), is determined by Fermi- 
Dirac quantum statistics because it reflects emergence 
from an environment where thermal equilibrium obtains. 
Moreover, the general relativistic analog to Liouville's 
Theorem assures us that as these neutrinos freely stream 
through spacetime, their proper phase space density re- 
mains unchanged and therefore continues to respect the 
Pauli principle. 

There is a second, peculiar effect arising from the fact 
that these relic neutrinos decouple in flavor eigenstates. 
A decoupled relic neutrino Ua can be regarded as a coher- 
ent superposition of mass states with common spacelike 
momentum p. As the universe expands and this momen- 
tum redshifts downward, there will come an epoch where 
the higher rest mass component of this flavor wavepacket 
will tend toward a nonrelativistic velocity, while the lower 
mass components speed along at nearly the speed of light. 
From this point onward the spatial scale of the coher- 
ent flavor wavepacket will grow rapidly with time. The 
net result is that at the present epoch, and at the red- 
shift z 1 epoch where neutrinos are falling into dark 
matter potential wells, the relic neutrinos are coherent 
structures with sizes considerably larger than the spa- 
tial scale of the gravitational potential wells of galaxies 
and, in some cases, comparable to the causal horizon 
size. Using the normal neutrino mass hierarchy scheme 
and mi = 1 meV, Fig.[4]shows how the size of a relic neu- 
trino flavor wavepacket would grow with the expansion 
of the universe. 

The more massive components of these neutrino fla- 
vor wave packets become nonrelativistic very early on, 
and it is interesting that the two different possibilities 
for the neutrino mass hierarchy yield quite different kine- 
matic histories for the relic neutrinos. For example, with 
mi = 1 meV in the normal neutrino mass hierarchy, the 
measured neutrino mass-squared differences imply that 
m2 ~ 8.7 meV, and ra-^ k, 49meV. As the wave packets 
for the three mass eigenstates redshift in the expanding 
universe, the components corresponding to each of the 
three mass eigenstates tend to become non-relativistic at 
three distinct epochs, z ~ 5, 50, and 300 respectively. 

On the other hand, the inverted neutrino mass hier- 
archy would yield a different history. The three neu- 
trino masses would be ma = 1 meV, mi « 49 meV, 
and m2 ~ 50 meV, which would imply that the light- 
est mass eigenstate would tend to be non-relativistic at 
z ^ 5 while the two heavier mass eigenstates would be- 
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FIG. 4: The separation of mass eigenstates as a function of 
redshift. Solid curve for ui—i/a, long dashed curve for 1^2 — fa, 
and medium dashed curve for ui — V2. The causal horizon is 
the short dashed curve. 

come non-relativistic at z ~ 300. 

If the lightest neutrino mass eigenvalue is <^ 1 meV, 
then the relic neutrinos might be coherent flavor eigen- 
states consisting of superpositions of relativistic and non- 
relativistic components at the epoch z « 2 vifhen these 
particles begin to fall into the gravitational potential 
wells associated with galaxies. Though the role of space- 
time curvature in quantum state reduction is not settled 
[nl, [13, [13] , an obvious set of questions is posed. 

First, does the process of "capturing" a neutrino into 
a gravitational potential well lead to flavor wavepacket 
de-coherence? Given the disparity between the rela- 
tively small local spacetime curvature scale and the flavor 
wavepacket size in this example, it is plausible that the 
gravitational tidal stress induces de-coherence. In this 
picture, the neutrino clustering/ capture process would be 



tantamount to a "measurement," with capture occurring 
when flavor wave function collapse yields a nonrelativistic 
mass-energy eigenstate. 

However, this line of reasoning calls into question the 
premise of coherent flavor eigenstates at the epoch when 
neutrinos are captured into, e.g., clusters of galaxies. Fig. 
|4] shows that the relic neutrino flavor wavepackets can 
have spatial extents that are comparable to the space- 
time curvature scale of the universe itself. In this pic- 
ture, the complete de-coherence history of the relic neu- 
trinos could be obtained only through solution of a fully 
covariant Dirac-like field development equation with all 
matter and dark matter distributions plus a prescription 
for de-coherence (wave function collapse). An alternative 
possibility is that de-coherence would never happen. 

In either case, unitarity implies that the comoving 
number density of neutrinos at a given mass is fixed. 
As a result, unless there is new physics associated with 
curvature-induced decoherence, the accretion history of 
neutrinos in potential wells should depend only on the 
neutrino masses and, in particular, the mass hierarchy. 

The on-going revolutions in observational cosmology 
and experimental neutrino physics overlap when it comes 
to the relic neutrinos, among the most abundant parti- 
cles in the universe. The fruits of this symbiotic over- 
lap include arguably the best probe of neutrino masses. 
The work presented here directly addresses these neu- 
trino mass issues, but also presents new insights into the 
role of quantum coherence and decoherence in the history 
of these relic particles. 
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